hep-th/9807109, IASSNS-HEP-98-65 



Theta Dependence in The Large N Limit 
Of Four-Dimensional Gauge Theories 



Edward Witten 

School of Natural Sciences, Institute for Advanced Study 
Olden Lane, Princeton, NJ 08540, USA 

The 9 dependence of pure gauge theories in four dimensions can be studied using a duality 
of large N gauge theories with string theory on a certain spacetime. Via this duality, one 
can argue that for every 9, there are infinitely many vacua that are stable in the large 
N limit. The true vacuum, found by minimizing the energy in this family, is a smooth 
function of 9 except at 9 = n where it jumps. This jump is associated with spontaneous 
breaking of CP symmetry. Domain walls separating adjacent vacua are described in terms 
of wrapped sixbranes. 
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1. The Problem 

In weak coupling, the dependence of four-dimensional gauge theories on the theta angle 
is computed via instantons. An instanton or anti-instanton contribution is proportional to 
exp (— 8iv 2 / g 2 ) exp(±i#). For example, in spontaneously broken gauge theories, in which 
instantons have a characteristic maximum size, and a characteristic effective coupling, the 
9 dependence is determined by an instanton expansion. 

For unbroken asymptotically free gauge theories, the situation is rather different. In 
such theories, at the classical level, instantons come in all sizes. The infrared behavior of 
the instanton gas is difficult to understand, and it is not clear what effective description 
should be used at long wavelengths to describe the theta dependence, or other aspects of 
the physics. 

One sharp way to pose the question is to consider the large iV limit of an SU(N) gauge 
theory ffl. The large N limit is an important avenue for understanding the dynamics of 
pure gauge theory, or gauge theory with a small number of light quark flavors, in four 
dimensions. The large iV limit is attained by taking N — > oo with A = g 2 N fixed. Thus 
the amplitude for an instanton or anti-instanton of definite size is weighted by a factor 
of exp(— 87r 2 iV"/A), and it appears that instanton effects would vanish exponentially for 
N — ► oo. However, there are a variety of reasons to believe that, because of infrared 
divergences, this is not so.0 For example, if light quarks are included, then the 9 dependence 
can be seen in current algebra HH, by reinterpreting some old pre-QCD computations 
0. One can show that if chiral symmetry breaking survives in the large iV limit, then so 
does the theta dependence. Moreover, the most plausible interpretation of how the solution 
of the U(l) problem of QCD fits into the 1/N expansion implies that in the pure gauge 
theory, the 9 dependence of the ground state energy is present to leading order! in 1/N fTj. 
These arguments suggest that, as in some two-dimensional models where somewhat similar 
questions can be asked ]8||]], the 9 dependence of pure gauge theory in four dimensions 
(with or without a small number of matter fields in the fundamental representation of the 
gauge group) is present in the leading order of the 1/N expansion. 

1 In contrast, A/" = 4 super Yang-Mills theory, which is scale-invariant rather than asymptoti- 
cally free, does not have these infrared divergences, and does have exponentially small ^-dependent 
effects, which can be computed via instantons Q. These instantons are related in the AdS/CFT 
correspondence to — 1-branes ||. 

2 An explanation of exactly what is meant by "leading order" will emerge below. 
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If so, one can draw an interesting deduction about the form of the theta-dependence 
|J. In any theory with N x N matrix fields $i, the large N limit is obtained by taking a 
Lagrangian of the form 

L($i) = NS(Qi-,w a ), (1.1) 

where S is independent of N and the w a are parameters such as bare masses and coupling 
constants. With the normalization in ( p. . 1|) , the large iV limit is obtained by keeping w a 
fixed as N — ► oo. 

Now instead, the most general renormalizable Lagrangian for gauge fields in four 
dimensions takes the form 

M 9 
L = —TrF^F^ + —e^TUT F a/J . (1.2) 

The normalization is chosen so that 9 is an angular variable. The general recipe of fll.lj ) 
would tell us to set 9 = Nip and keep ip fixed for N oo. In the large N limit, the vacuum 
energy E is proportional to iV 2 (as the number of degrees of freedom is of that order). So 
we expect E(?p) = N 2 h{ip) for some function h which should have a limit as N — ► oo. In 
terms of 0, that means 

E{9) = N 2 h(9/N). (1.3) 



In addition, E must obey 



E{6) = E(9 + 2n). (1.4) 



These conditions are, however, practically incompatible: a smooth function of 9/N cannot 
be invariant to 9 —>■ 9 + 2n unless it is constant. 

The most plausible way out seems to be || that E{9) is a multibranched function 
because of many candidate vacuum states that all become stable (but not degenerate) for 
N = oo. Such behavior occurs in many two-dimensional models |H| (including some 
with 1/N expansions that raise issues like those we are discussing here). In the k th vacuum, 
the energy would be 

E k (9) = N 2 h((9 + 27tk)/N). (1.5) 

The truly stable vacuum would be found, for each 9, by minimizing E^ with respect to k. 
The actual vacuum energy would be therefore 

E(9) = N 2 mmh((9 + 2nk)/N). (1.6) 

k 



This function is periodic in 9, but (if h is not constant) it is not smooth - at some value 
of 9 there is a jump between two different "branches." 

Under a CP transformation, one has 9 — ► —9. So in particular CP is a symmetry if 
and only if 9 equals or tv. Hence h{9) = h(—9). Note that CP acts by k — > — at 9 = 0, 
and by /c — > — 1 — k at # = it. 

Moreover, E{9) has its absolute minimum at 9 = 0, because precisely at 9 = the 
integrand of the Euclidean space path integral is real and positive.! If the vacuum is 
unique at 9 = 0, then the minimum in (|1.6| ) occurs for k = (otherwise k and —k would 
both contribute). Moreover, one expects that d 2 h/d9 2 ^ at 9 = because of arguments 
involving the U(l) problem in the theory with quarks 0, or simply because of the absence 
of a symmetry that would make this quantity vanish. If so, we can set h(9) — C9 2 + . . ., 
where C is positive and the higher order terms do not contribute to ( |1.6| ) in leading order 
in 1/N. Thus, one conjectures that the large iV structure of the vacuum energy is 

E(9) =Cmm(9 + 2nkf + 0(1/N), (1.7) 

k 

with some constant C. This function exhibits a nonanalyticity at 9 — it, which we associate 
with a jump between two vacua (with k = and k = —1) and the spontaneous breaking 
of CP invariance. 

2. The Computation 

In what dynamical approximation to gauge theory can one hope to check the ideas 
that were just reviewed? All known approaches to the dynamics of four-dimensional gauge 
theory in which one can exhibit any of the difficult properties like confinement and the 
mass gap involve replacing the theory by a simpler theory which is hoped to be in the same 
universality class. There are several candidates for what the simpler theory can be. One 
candidate is lattice gauge theory. This framework makes it possible to compute a great 
deal, especially upon using computer simulation to get away from the strong coupling 
limit. But it is not very convenient for discussing the 9 dependence, particularly if one 

3 At 9 = all contributions to the path integral receive positive weights; at 6 ^ 0, this is not 
so because the instanton factor is e %e . The Euclidean space path integral in volume V computes 
exp(— VE(6)), so E(Q) is minimized by maximizing the Euclidean space path maximal, which 
happens when the weights are all positive. 
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wishes to probe issues in which the invariance under 9 — > 9 + 2n is important. Another 
possibility is to consider the realization of four-dimensional gauge theory via M-theory 
fivebranes ||11|| , where again, at the cost of replacing the theory of interest by a simplified 
version that is hopefully in the same universality class, one can demonstrate the mass gap 
and confinement. We will instead study the problem in yet a third framework in which 
one can exhibit the mass gap and confinement in the context of a simplified version of 
four-dimensional gauge theory. This involves a circle of ideas connected with the 

correspondence between conformal field theory and quantum gravity on Anti-de Sitter 
space 

To get to the specific issues of interest here as quickly as possible, we will begin as in 
section 4 of [15| with Type IIA superstring theory on M = R 4 x S 1 x R 5 , with N parallel 
fourbranes whose worldvolume is V = R 4 x S 1 x x; here a; is a point in R 5 . We assume that 
the "spin structure" is such that fermions change sign in going around the S 1 . Then the 
theory on the branes is at low energies a pure U(N) gauge theory in four dimensions. It 
follows from the general AdS/CFT correspondence that the large N behavior of the SU(N) 
part of this gauge theory can be studied by studying weakly coupled string theory on the 
supergravity solution X which these branes generate. Topologically X = R 4 x D x S 4 , 
where D is a two-dimensional disc. The change in topology from M to X is crucial (along 
with the fact that A is a smooth manifold without branes), both in the explanations of 
confinement and the mass gap in Jl5[ and in the discussion below of the 9 dependence. 
The metric of X is 

ds 2 = f ^A 3 |>T + J^Att (a 2 - 1) + f ^A^V + fvXdnl (2.1) 

Here x % are coordinates on R 4 , A and ijj (with 1 < A < oo, < ijj < 2n) are polar 
coordinates on D (note that A = 1 is the origin of the polar coordinates, the "center" of 
D) , and dO| is the metric of a unit four-sphere. 77 is a parameter which determines how far 
one is from conventional four-dimensional gauge theory; for r\ » 1, the string theory on 
X can be studied via long wavelength supergravity, while asymptotically-free gauge theory 
is expected to emerge in the opposite limit 77 — > 0. In the present paper, we work in an 
approximation of long wavelength supergravity, so to compare with gauge theory, we must 
assume that the system has no phase transition as a function of r\. 

How can we include 9 in the formulation of gauge theory via fourbranes? This can be 
done quite simply by including the U{1) gauge field that arises in Type IIA superstring 
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theory from the Ramond-Ramond sector. We will in this paper denote that field as a, 
and its field strength as = didj — djdi. Let us reconsider Type IIA superstring theory 
on M = R 4 x S 1 x R 5 with the wrapped fourbranes of worldvolume V. The low energy 
worldvolume effective Lagrangian of the fourbranes has a term 

f TrF AF . . 

(The most familiar manifestation of this term is that instantons on the fourbrane are 
charged with respect to a - they carry zerobrane charge.) Here F is the U(N) field 
strength. We now modify the Type IIA vacuum so that / = 0, but 



a = 9 a (2.3) 



is possibly non-zero. (The left hand side is gauge-invariant modulo 27rZ, so we interpret 
6 a as an angle.) At low energies in four dimensions, ( |2.2| ) reduces to a theta term in the 
gauge theory action, and the four-dimensional Yang-Mills theta angle is 

e = e a . (2.4) 

What we have done so far is just to learn how to include 9 in the fourbrane description 
of four-dimensional gauge theory. Now, we go over to the dual description by supergravity 
(or string theory) on X. In doing so, we must bear in mind that the parameters of the 
theory are determined by specifying the Type IIA vacuum far away from the branes, that 
is at large A, and then the behavior at small A is determined by the supergravity equations 
on X (or, if r\ is small, the full string theory on X) and encodes the behavior of the gauge 
theory. For example, in the original description with branes on M we assumed that / = 0. 
In the dual description on X , the analogous statement is that / = for A — > oo. Likewise, 
( |2.3|) should be interpreted to mean that J gl a = 9 a = 9 if the integral is taken at large 
A. If we combine these conditions (plus Stokes' theorem f D f = lim^oo J s i we learn 
that j D f = 9 a = 9 mod 2ttZ. The 27rZ indeterminacy arises because the left hand side is 
a well-defined real number, but 9 is an angle. Hence 

f = 9 + 2nk (2.5) 



D 



for some integer k. 
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J D f is of course defined as J D dXdipfx^. 
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Maxwell's equations for the / field have a normalizable zero mode in which the only 
nonzero component is 

= ^f- (2-6) 
The normalization has been chosen so that f D f = 2n. Hence, we can find a solution of 



(2.5) obeying Maxwell's equations in the simple form 



Av, = (0 + 27rfc)4r (2.7) 

TV A' 

The back reaction on the geometry produced by this / field is negligible in the limit that 
N goes to infinity with fixed 9 + 2nk. The reason for this is that the classical action of the 
spacetime is of order N 2 (like the vacuum energy of the large N gauge theory to which it 
is dual). The kinetic energy of the / field is of relative order 1/N 2 as it is simply 



/ 



with no factors of N or of the string coupling constant. So to lowest order in 1/N, we obey 
Q2.5Q and the classical equations of motion simply by solving for / in the fixed spacetime 
X. 

The 9 and /c-dependent part of the vacuum energy is found by evaluating ( |2.8|) with 
/ as given in ( |2.7|) . It thus takes the form foreseen in the introduction: 

E k (9) = C(9 + 2nk) 2 , (2.9) 

with some positive constant C that is independent of N. The vacuum energy for given 9 
is obtained by minimizing this with respect to k, 

E(9)=Cmm(9 + 2nk) 2 . (2.10) 

k 

We have obtained precisely the structure anticipated in the introduction. For given 9, 
there are infinitely many vacua, labeled by the choice of an integer k. The vacuum energy 
E{9) is a smooth function of 9 except at 9 = it, where there is a jumping between two 
solutions (k = and k = — 1). This jumping represents spontaneous breaking of CP at 

9 = 7T. 

Finally, we should check that the vacua labeled by k are all stable for N — > oo, and find 
a framework for estimating their lifetime for finite N. The essential point is to describe the 
domain wall separating vacua with adjacent values of k. The decay of a k- vacuum involves 
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nucleation of a "bubble" with a smaller value of \6 + 2irk\ and hence a lower energy density; 
this bubble is bounded by a domain wall. If the energy per unit area of such a domain 
wall is large for iV — > oo, then the k- vacua have lifetimes that go to infinity for N — > oo. 

In fact, the domain wall is constructed simply by compactifying a Type IIA sixbrane 
on the S 4 factor in X = R 4 x D x S 4 . Let w be a point in D and let C be a codimension 
one surface in R 4 given by, say, x 3 = (with x 3 one of the space coordinates) . Consider a 
sixbrane whose worldvolume is Q = C x w x S 4 . The value of k jumps by ±1 in crossing 
such a sixbrane (the sign depends on orientations and on the direction of crossing the 
sixbrane) . The defining property of a Type IIA sixbrane is that if E is a two-surface that 
has linking number 1 with the sixbrane worldvolume, then f E f = 2n. In the present case, 
we can take E = D\ — D2, where D\ and D2 are copies of D that are respectively to the 
"left" or "right" of C.i So 

/ /- / / = ±2tt, (2.11) 

J D 1 J D 2 

and this means that k jumps by ±1 in crossing the sixbrane. The domain wall is thus not an 
ordinary soliton, as one might naively have thought, but a .D-brane. In particular, color flux 
tubes associated with quark confinement can terminate on the boundaries between different 
vacua, just as they can terminate |nj on chiral domain walls in M = 1 supersymmetric 
gauge theory in four dimensions. (Both kinds of domain wall involve a 2tt jump in 9.) 

Now, the energy density of a Type IIA sixbrane is for weak coupling of order l/\ s t 
(X s t is the Type IIA string coupling constant), and as a result, in the large N limit it is 
of order N. An instanton describing the decay of a /c-vacuum can be constructed as a sort 
of sixbrane bubble and has an action that grows as a power of N. So the lifetime of a 
k- vacuum, for any given k, is exponentially long for iV — > 00. 

This work was supported in part by NSF Grant PHY-9513835. I would like to thank 
S.-J. Rey for comments. 



5 In other words, Di = yi x D x z, where z is a point in S 4 , and the yt are points in R 4 to the 
left or right of C. 
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